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One of the most interesting development trends of a modern cosmology is the analysis of models 
of a modified gravitation. Without exaggeration it is possible to say that Sergei Odintsov is one 
of the leaders of this direction of researches (see [lj]). This article is dedicated to cosmologies with 
variable speed of light (VSL) - models, which one can consider as a particular case of models of a 
modified gravitation. 

In quantum cosmology the closed universe can spontaneously nucleate out of the state with no 
classical space and time. The semiclassical tunneling nucleation probability can be estimated as 
P ~ exp(— cv 2 /A) where a=const and A is the cosmological constant. 

In classical cosmology with varying speed of light c(t) it is possible to solve the horizon prob- 
lem, the flatness problem and the A-problem if c = sa n with s=const and n < —2. We show 
that in VSL quantum cosmology with n < —2 the semiclassical tunneling nucleation probability is 
P~ exp(— (3 2 A k ) with /3=const and k > 0. Thus, the semiclassical tunneling nucleation probability 
in VSL quantum cosmology is very different from that in quantum cosmology with c=const. In 
particular, it can be strongly suppressed for large values of A. In addition, we propose two instan- 
tons that describe the nucleation of closed universes in VSL models. These solutions are akin to 
the Hawking- Turok instanton in sense of 0(4) invariance but, unlike to it, are both non-singular. 
Moreover, using those solutions we can obtain the probability of nucleation which is suppressed for 
large value of A too. We also discuss some unusual properties of models with n > 0. 

PACS numbers: 98.80.Cq;98.80.-k 



I. INTRODUCTION 

One of the major requests concerning the quantum 
cosmology is a reasonable specification of initial condi- 
tions in early universe, that is in close vicinity of the 
Big Bang. There are known the three common ways to 
describe quantum cosmology: the Hartle-Hawking wave 
function the Linde wave function and the tun- 
neling wave function Q. In the last case the universe 
can tunnel through the potential barrier to the regime 
of unbounded expansion. Following Vilenkin [j| lets con- 
sider the closed (k = +1) universe filled with radiation 
(w = 1/3) and A-term (w = —1). One of the Einstein's 
equations can be written as a law of a conservation of the 
(mechanical) energy: P 2 + U (a) — E, where P = —act, 
a(t) is the scale factor, the "energy" E — const and the 
potential 




FIG. 1: Vilenkin potential with c =const. 
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where c is the speed of light; see Fig.l The maximum 
of the potential U(a) is located at a e = y / 3/2A where 
U{a e ) = 3c 2 /(4A). The tunneling probability in WKB 
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approximation can be estimated as 
t 2c 2 i" 



P ~ exp 
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da^U(a) - E , (1) 



where < Oj are two turning points. The universe can 
start from a = singularity, expand to a maximum ra- 
dius a! i and then tunnel through the potential barrier to 
the regime of unbounded expansion with the semiclassi- 
cal tunneling probability |T]). Choosing E = one gets 



= and 



/3/A. The integral in fl} can be calcu- 



lated. The result can be written as 
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For the probability to be of reasonable value, for example 
P= 1/e ~ 0.368, one has to put A - 0.3 x 10 65 cm" 2 (see 
©). In other words, the A-term must be large. How- 
ever, despite this problem, we does acquire one prise: 
Once nucleated, the universe immediately begins a de 
Sitter inflationary expansion. Therefore the tunneling 
wave function results in inflation. And the A-term prob- 
lem, which arises in this approach is usually being gotten 
rid of via the anthropic principle. 

Now let's return to the Fig. 1. We have two Lorentzian 
regions (0 < a < a' it a > ai) and one Euclidean region 
{a[ < a < ai). The second turning point a = ai cor- 
responds to the beginning of our universe. If A = 
then U(a) has the form of parabola and we get only one 
Lorentzian region. In this case, the universe can start 
at a = 0, expand to a maximum radius and then rec- 
ollapse. If E — > then the single Lorentzian region is 
contract to the point. This, of course, comes to an agree- 
ment with the tunnelling nucleation probability: P — > 
as A — > 0. In this article, however, we'll show that in 
quantum cosmological VSL models the situation can be 
opposite, viz: the probability to find the finite universe 
short after it's tunneling through the potential barrier is 
P~ exp(-/3(n)A a(n )) with a(n) > and 0(n) > when 
n < — 2 or for — 1 < n < — 2/3. After the tunneling one 
get the finite universe with "initial" value of scale fac- 
tor dj ~ A -1 / 2 , so the probability to find the universe 
with large value of A and small value of a* is strongly 
suppressed. The reason of this lies in the behavior of 
potential U(a), which, for the case A — * 0, is being trans- 
formed into the hyperbola, located under the abscissa 
axis. As a result, such a universe can start at a ~ 
the regime of unbounded expansion. Therefore, we get 
the single Lorentzian region which is not contract to the 
point at E — » 0. 

This new property of VSL quantum cosmology will be 
discussed in the Sec. II for the case w = 1/3. But there 
arouse the two new questions which have to be answered. 
First at all, the effective potential in VSL models can be 
unbounded from below at a — > 0. What possible mean- 
ing of such potentials can be? The second question is 
the geometric interpretation of the quantum creation of 
a Universe with varying speed of light. We know that 



universe can be spontaneously created from nothing (in 
model where c = const) and this process can be described 
with the aid of the instantons solutions possessing 0(5) 
(if V{4>) has a stationary point at some nonzero value 
4> — 4>o = const) or 0(4) (as Hawking- Turok instanton 
6]) invariance. So, what can be said about instantons in 
the VSL models? 

The whole plan of the paper looks as follows: in 
the next Section we'll consider the simplest VSL model: 
model of Albrecht-Magueijo-Barrow. Then we show that 
in framework of tunneling approach to quantum cos- 
mology with VSL the semiclassical tunneling nucleation 
probability can be estimated as P ~ exp(— f3 2 A k ) with 
/3=const and k > 0. All corresponding calculations will 
be done for the case of the universe filled with radia- 
tion (w = 1/3) and vacuum energy. The case of any 
w will be considered in the Sec. Ill, where we will dis- 
cuss the problem of potentials unbounded from below at 
a — > 0. We'll use the following naive procedure: By anal- 
ogy with Schrodinger equation we apply the Heisenberg 
uncertainty relation to the Wheeler-DeWitt equation in 
order to find the potentials which admit the ground state, 
even being unbounded from below at a — ► 0. 

In the Section IV we'll propose the non-singular in- 
stanton solutions possessing only 0(4) invariance (so the 
Euclidean region is a deformed four sphere) . These solu- 
tions can in fact lead to inflation after the analytic contin- 
uation into the Lorentzian region. We will discuss these 
results in Sec. V. Some unusual properties of VSL mod- 
els with n > (including "big rip" and "big trip") are 
discussed in Appendix. 



II. ALBRECHT-MAGUEIJO-BARROW VSL 
MODEL 

Lets start with the Friedmann and Raychaudhuri sys- 
tem of equations with k = +1 (we assume the G=const): 
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where a = a(t) is the expansion scale factor of the Fried- 
mann metric, p is the fluid pressure, p is the fluid density, 
k is the curvature parameter (we put k = +1), A is the 
cosmological constant, cq is some fixed value of speed of 
light which corresponds some fixed value of scale factor 
a . 

Remark 1. In [7] George F.R. Ellis and Jean-Philippe 
Uzan has pointed out at the fact, that the system of 
Friedmann- Raychaudhuri equations (system as de- 
noted in our article) is not consistent from the point 
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of view of a field theory. This is because the system 
© couldn't be derived via the standard variation pro- 
cedure, despite the contrary claims, given in [7]- [10] (see 
the references in 0). Variation of action (36) (see 0) 
leads instead to the system of equations (39)- (42) (see 
0), and not to the Friedmann-Raychaudhuri. There- 
fore, since our article is initially based on ^ it seems 
that the present article is just a formal exercice based 
on unreliable grounds (Eqs. (J3])) so that its conclusions 
cannot tell us much about the effect of varying constants 
in Quantum cosmology. 

We don't think so. The fact, that the equations (39)- 
(42) (of 0|) are distinct from set ([3]) doesn't immediately 
mean that the system (39)-(42) (of 0) is correct while 
<j3j> is not. Moreover, even assuming that George F.R. 
Ellis and Jean-Philippe's claim is indeed true and it is 
impossible to get fl3) from (36) (see 0), this would not 
have any consequences for our further results, barring the 
instanton chapter. In our article (except the mentioned 
chapter) we use neither the action (36) (of [7|) nor the 
fact that §5§ can be derived from (36). Contrary to that, 
we just use © as a basic phenomenological model. This, 
of course, can be regarded as a flaw of a model, but same 
can be said about the system (39)-(42) (of [7j]), which 
is also but a phenomenological model, based on assump- 
tion, that the VSL models are in essence the particular 
example of the scalar-tensor theories (see for example 

S). 

Hence, in order to be able to make the preference of 
(39)- (42) (of 0) over equation §5§ (or vice verse) we have 
to gain the better comprehension of the physical princi- 
ples, lying behind the variability of the speed of light (in 
assumption, that it really changes at all, which is still 
far from obvious!..). That, in turn, would be possible 
only after we will understand the origins of " c" . It has 
been shown in reviewing part of 0, that it is possible 
to define not just one, but a couple of "speeds of light" 
(cst, Cem, ce). Beyond any doubts, there should exist 
a fundamental physical reason for those (generally differ- 
ent) values being perfectly equal. Such reason can only 
be discovered in the future fundamental physical theory 
(string theory?). Only such theory would verify which 
equations are true: (39)-(42) (in 0), ([3]) or maybe some 
other ones. 

Returning to system (39)-(42) (in 0), which may be 
considered as a more reliable one than Friedmann and 
Raychaudhuri system ([3]) , we shall in turn remind about 
some of the assumptions, that has been used by George 
F.R. Ellis and Jean-Philippe's in order to get to (39)-(42). 
For example, authors consider the value of ce (which is 
used in Einstein's equation in a form 8ttG/c 4 ) as a vari- 
able, while assuming that the value of cst , that appears 
in the integral 

ds 2 = -c 2 ST dt 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 2 ) 2 , 

satisfies c$t = 1 (see page. 8, in the last line above 
the formula (39) in 0). But if we'll assume that all 
"speeds of light" are in fact the representations of just 



one value (an assumption, that we strongly believe to 
be true), then the value of c = i/? 1 / 4 will require to be 
multiplied by the lagrangian (d A x — cdtd 3 x), including 
the various fields of matter and shall also be taken into 
account in the multiplier ^fg. With all this in mind, it 
becomes really doubtful that the VSL model (36) ([7j) 
will remain being equal to some scalar-field theory. On 
the other hand, if we'll insist on difference between the 
values of c$t and ce, while requiring the variability of 
both of them, it doesn't seem impossible to choose the 
cst in a way, that leads to results, similar to the ones, 
received in [7]- [10] (see the reference list of 0), which is 
sufficiently what we need for the IV chapter of our work, 
devoted to instantons (and is necessary for this chapter 
only!). Of course, this might look quite artificial, but it 
just serves to show that both (39)- (42) of [7| and ^ shall 
be regarded as phenomenological models. 

Therefore, it appears, that studying of equations ([3|) 
are no more (or less) well-grounded, then studying of 
equations (39)- (42) of [7|. As follows from our work, sys- 
tem ([3]) can be quite interesting in the quantum cosmol- 
ogy, since it allows for an unexpected solution of a cosmo- 
logical constant problem. This, with regard to the overall 
difficulty of the problem, can be considered as an addi- 
tional indication that the model ([3]) does indeed deserve 
the further examinations. We aren't making any further 
claims; it is of course quite possible that the system (39)- 
(42) of 7], if being applied to the quantum cosmology is 
also capable of producing of the interesting results, but 
that is a topic of a different research. 

In conclusion, let us note that the equations (39)- (42) 
of [7j are in fact completely incompatible with our model. 
Our article is based on assumption that there exist a 
one-to-one correlation between c and the scalar factor a: 
c = sa n . In other words, the speed light (or ifj = c 4 
from 0) is a function of a Ricci scalar R, and, hence, 
is NOT an independent variable. But then we come to 
contradiction with the basic assumption of , which one 
basically allowed to get the system (39)- (42). Therefore, 
it would even formally be a mistake to use the (39)-(42) 
of in a way, similar to the one, that we adopted for 
([3]) - at least without special pre-given restrictions. 

Remark 2. One can consider more general model. 
Let 
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with n > 0, m > 0, r > 0, s > 0. In this case the flatness 
and A-term problems can be solved if 



n > to, n > 1, m > 



3(w + l) 
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for the weak energy condition is valid: w > — 1. If a — > oo 
then the Hubble root 
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so we have the big rip singularity without any phantom the expression p which has the form: 
fields (see Appendix A). All this can be obtained from 
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After that extensive digression let's return to the equa- 
tions ([3]). Using these ones 
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P = [P + ~l ^ a r> 2 — • 

Choosing w — 1/3 one can solve ((5]) to receive 
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where M > is a constant characterizing the amount of 
radiation. It is clear from the © that the flatness prob- 
lem can be solved in a radiation-dominated early universe 
by an interval of VSL evolution if n < — 1, whereas the 
problem of A-term can be solved only if n < —2. The evo- 
lution equation for the scale factor a (the second equation 
in system ([3])) can be written as 



p 2 + U(a) = E, 



(7) 



where p = —ad is the momentum conjugate to a, E = 
8?rGM/3 and 
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The potential ([5]) has one maximum at a = a, 
v/3/(2 A) such that 

U e = U{a e ) 



s 2 3 „+l 



2«+ 1 A»+ 1 (n+l)(n + 2) : 



(9) 



so U e > if (i) n < —2 or (ii) n > —1. The first case 
allows us to solve the flatness and "Lambda" problems. 
The surplus dividend of the model is the presence of fi- 
nite time region during which universe has accelerated 
expansion. 



The equation ([7]) is quite similar to equation for the 
particle of energy E that is moving in potential (|10p . 
hence the universe in quantum cosmology can start at 
a ~ 0, expand to the maximum radius a[ and then tunnel 
through the potential barrier to the regime of unbounded 
expansion with "initial" value a = dj. The semiclassical 
tunneling probability can be estimated as 



(11) 
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P~ exp -2 / | p(a) \da , 



c 2 (t) 



8ttGH 



\p(a)\, \ P (a) |= y/U m (a)-E, 



where E < U e . It is convenient to write E — U e sm 2 9, 
with < 6 < n/2. 

For the case E <g; U e one can choose 
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and evaluate the integral (fTTj) as 
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where 
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A. The semiclassical tunneling probability in VSL 
models with n < —2: the case E -C t/ e 

One can choose n = —2 — m with m > 0. Such a 
substitution gives us the potential ([8} in the form 
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The integral ([14]) can be calculated for the m G Z. For 
example 
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(3 + cos 2 6) (191 - 78 cos 2 9 + 15 cos 4 6) \/6 + 2 cos 2 ~ 0.148 + O(0 6 ), 



r 



h ~ 0.025, 



h ~ 0.007, J 4 ~ 0.002, 



and so on. One can further show that I m (9) > at 
< 9 <C 1. Thus, it is easy to see from 1)13)1 that the 
semiclassical tunneling probability P — > for large values 
of A > and P -> 1 at A -> 0. 

Note, that the case c=const can be obtained by sub- 
stitution m = — 2 into the (fTB")) . Not surprisingly, this 
case will get us the well known result P ~ exp(— 1/A) 
(see @). 



B. The semiclassical tunneling probability with 

n < —2 and n > — 1 



z is dimensionless quantity and z[, Zi are the turning 
points, i.e. two real positive solutions of the equation 
F m (z,9) — for the given 9 (it is easy to see that the 
equation F m {z,9) — does have two such solutions at 
< 9 < tt/2). 

If m is the natural number then the expression (|15|) 
has a more simple form. For example 



Pi ~ exp [ ^ I — yj (z 2 - z[ )(zf - z 2 ) , 
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with 



In the case of general position the semiclassical tunnel- 
ing probability with n = — 2 — m has the form 



P m ~ exp 



where 
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This expression can be calculated exactly: 
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where ^ 2 = cos6>/cos 2 (6»/2), A = cot(6>/2), n and K are 
complete elliptic integral of the first and the third kinds 
correspondingly [§]. 

Similarly, P ~ exp(— S 1 ), with 
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regime of unbounded expansion is strongly suppressed for 
large values of A and small values of the initial scale factor 
at = v / 3~/(2sin(6 , /2)\/A). In other words, overwhelming 
majority of universes which are nascent via quantum tun- 
neling through the potential barrier ([8]) have large initial 
scale factor and small value of A. 

Now, let us consider the case (ii), when n > — 1. The 
"quantum potential" has the form 



U(a) 
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and so on. 

Therefore the probability to obtain (via quantum tun- 
neling through the potential barrier) the universe in the 



where m = n + 1 > 0. The points of intersection with the 
abscissa axis a are a = and ai — y/3(m + 1) /2Am. 
Choosing E — in equation ^} and substituting (fT7| 
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into the (jlip we get 

/ s 3 A (i- 3m) /2 n \ 

p ~ exp (- 4^ A z p-s^Ti) dz )' 

with zi = ^/3(m + l)/2m (The starting value 2 = 
means that the Universe tunneled from "nothing" to a 
closed universe of a finite radius a\ = Z\ /VI.). Thus, we 
have the same effect asif0<m<l/3. 

C. Peculiar cases with n = — 1 and n = —2 



Therefore, U e > if a < 2 e - 3 / 4 /\/6 - 0.386. Choosing 
a = 0.286 and E <C C/ e gets us the turning points z- ~ 
0.77 and z t ~ 2.391. 

At last, the semiclassical tunneling nucleation proba- 
bility is 

/ s 3 a 2 r zi dz r~\ 4 ~ ~ r\ 

/ 0.084s 3 A 2 \ 



At last, lets consider the cases of n — — 1 and n = — 2. 
The formula (fT5|) is not valid in these cases (m = — 1 and 
to = 0) so we shall consider these models separately. 

If n = —1 (to = —1) then 
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where a* is constant and [a*]=cm. The potential (|18|) 
has one maximum at a — a e = y // 3/(2A) such that U e = 
U(a e ) = 2s 2 log(a e /a*), so if a e > a* then £/ e > 0. We 
choose a* = A^ 1 / 2 . This gives us J7 e = 0.41s 2 > 0. For 
the case E <§; U e the semiclassical tunneling nucleation 
probability is 
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where the turning points are z[ = 0.721, Z{ = 1.812. As 
we can see from the (|19|) . when n = — 1 we receive the 
aforementioned effect again. 
If n = -2 (m = 0) then 
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We choose a* = 1/ (a-\/A) , where a is a dimensionless 
quantity. Thus 



U(a) 



1 4A 



log (aaVlCj 



(20) 



The maximum of potential (|2"P|) is located at the same 
point a e and 



s 2 A / , /9a 4 
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III. AN EXISTENCE OF GROUND STATES 
FOR SINGULAR POTENTIALS 

It is interesting to ask: what can be said about the 
value of U(a) when a = 0? What can be the possible 
meaning of the potential which at a — > is unbounded 
from below? It seems that such universe is able to just 
roll down towards small values of a (where the potential 
is tending to minus infinity) instead of any tunneling to 
large values. 

This situation can in fact be alleviated if the consid- 
ered potential U(a) has the ground state. Indeed, one 
can imagine the fictitious particle with some energy and 
coordinate a(t) in the potentials © rolling down to- 
wards small values of a. The main problem is: whether 
the quantum mechanical energy spectrum of U(a) is un- 
bounded below? If not, then it does admit the ground 
state and hence can have the physical meaning. 

To find such a potential lets suppose that our fictitious 
particle is located in a small region a near the singularity 
a = 0, with the momentum P. We will consider the case 
of arbitrary w. In this case the flatness problem can be 
solved in early universe by an interval of VSL evolution 
if n < n fl {w) = —(1 + 3w;)/2, whereas the problem of 
A-term can be solved only if n < n A (w) — n fl (w) — 1 = 
-3(w + l)/2. 

One can use the Heisenberg uncertainty relation as 

Pa~ (87rGn) (1 -"/' (t0))/2 c (3 "/' (tu) - 1)/2 . (21) 
Using (JUl), and © (or 0) one get for the a -> 
E = P 2 + U(a) -> 
Z 2 s 2 n Jl (w) 



a 2 - n{3n f (tu) - x) (n fl (w) - n)a 2{n n (™>-»0 ' 

where Z 2 = {SkGK) 1 ^! 1 s 3n n~\ and 

n <n fl (w) < 0. (22) 
Therefore the energy spectrum will be bounded below if 



(3n + 2) (n fl (w) - l) < 0. 
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FIG. 2: The ground state exists for w and n from the interior 
of the triangle ABC. 



Here we get the following Euclidean equations: 
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where primes denote derivatives with respect to r. 

At the next step we represent the potential V in fac- 
torized form 



V = F(a)U((j>). 



(26) 



Indeed, lets for example consider the power-low potential 
~ <j) k . If the coupling A is dimensionless one then we get 



V 



-G k ' 2 - 2 c 7 - 2k 4> k . 
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and Ij22p are valid. This situation is represented graphi- 
cally on the Fig. 2. It is easy to see that the conditions 
of ground state existence can be satisfied if n A (w) < n < 
n fl (w). (The additional restriction w > —1/3 is just a 
condition of existence of maximum of the potential U (a). 
For more detailed examination of the general case see 
[To|). In the case of universe filed with radiation (as 
above) one get —2 < n < —1. 



Since c = sa n then in the simplest case we come to (|26|) . 

Let = 0o = const be solution of the (|25j) . (Note, that 
we don't require the 0o to be the extremum of potential.) 
Using the first equation of system (|2"5)) and (f2"6")) we get 
the equation for the F(a), 



dF(a) 
da 



2n 



F(a) 



3ns 



,4n.-3 



ns A 
47rGC/ n 



- 1 , (27) 



where Uq — U{4>o) — const. The integration of the 
results in 



IV. INSTANTONS 



If we are going to describe the quantum nucleation of 
universe we should find the instanton solutions, simply 
putted as a stationary points of the Euclidean action. 
The instantons give a dominant contribution to the Eu- 
clidean path integral, and that is the reason of our inter- 
est in them. 

First at all, lets consider the 0(4)-invariant Euclidean 
spacetime with the metric 



F(a) 



C 



3ns 



8nG(l - n)U 



,2(7.-1) 



8irGU 



where C is the constant of integration and by assump- 
tion n / -1 and n / 0. Substitution of (|28|) into the 
second equation of the system (|25]l transforms it into the 
the model of nonlinear oscillator the integration of which 
result in 



where 



— +u{a) = 0, 



(29) 



ds 2 = c 2 (r)dT 2 



! (r) (dtf 



(24) 



In the case c = const one can construct the simple 
instantons, which are the 0(5) invariant four-spheres. 
Then one can introduce the scalar field <fi, whose (con- 
stant) value 4> = 0o is chosen as the one providing the 
extremum of potential V ((/>). The scale factor will be 
a(r) = i? _1 sin Ht and after the analytic continuation 
into the Lorentzian region one will get the de Sitter space 
or inflation. Many other examples of non-singular and 
singular instantons were presented in (Tlj 

Now, lets consider the VSL model with scalar field. 



u(a) 



2(1 -ny 



(30) 



with ui 2 = 8ttGU C/(3s 2 ) and with the choice C > 
made. We can see that for c = const (i.e. n = 0) the 
equation (|30[) turns out to be just the usual harmonic 
oscillator and we come to the well-known 0(5) solution 
(but in this case 0o must be the stationary point of V). 

The equation (|2"9")l naturally describes the "movement 
of a classical particle" with zero-point energy in mechani- 
cal potential ([30)) . Depending on value of n this potential 
can take one of four distinct forms (excluding the well- 
known classical case n = 0, which lies beyond the scoop 
of this article). 
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FIG. 3: The potential u(a). Instanton exists in region < 
a < ai but only if n > —1/5. 



Case 1: n < 0. Potential u(a) has the form, depicted 
on the Fig. 3. Here we have one Euclidean (0 < a < ai) 
and one Lorentzian (a > ai) regions where 



«i 



l/(l-n) 



(31) 



On the bound between Euclidean and Lorentzian regions 
(a = ai) we have a' = 0. 

This mechanical potential is unbounded from below at 
a — > 0. With this in mind, we'll have to ascertain that 
the Euclidean action for our solution will stay finite. The 
gravitation action has the form 



d A x- 



8nG 



VgR- 



We are using the dimensionlcss variables x° — cqt/oo, 



ip and so on. Calculating R we get 



R 



(i 



2„2 



eg - (^) 2n ((l-n) fl ' 2 + aa") 



(32) 



so we do have the potential divergence at a = 0. Multi- 
plying (|32p on the and c 3 and using the equation of 
motion we get the expression: 



R^c 3 ~ 6c (2 - n.)w 



,2n+3 



m— I i 1r< 

i t ~ n a Q 



(33) 



where the most dangerous multiplier factor is a 1+4n . But 
if— l/4<n<0 then the Euclidean action becomes finite 
and therefore, we end up with the legitimate gravitation 
instanton. In a similar manner, using (|26p and (|28[) we 
get for the scalar field (in dimensionless x^): 

l-3n 



3nc%a 1+5ri 



(n — l)a 



2/i 



A^Q 3 + 5r ' 
3n 



therefore the instanton exists for n > —1/5. This demand 
is more powerful than what we got for the gravitation 
instanton where n > —1/4 (see (|33j)). 

Case 2. < n < 1. Here the potential w(a) suffers 
no singularity at a = 0, but w(0) = 0. Also this potential 
has a minimum at 



a 



!/(!-„) 



and is equal to zero at (|3ip . hence, once again creating 
one Euclidean and one Lorentzian regions, separated by 

©. 

Case 3. n = 1. This case is somehow special, since 
for such n the solution of (|27p shall be 



F{a) =a 2n [C 



3s 4 



4:7tGU 



In a 



s 4 A 
8ttGU ( 



instead of (|2"5| . and hence, the equation of ([50)1 shall be 
substituted by 



(34) 



u(a) = a ( — s In a 



It is easy to see that this function has two zeros (at a\ = 
and a,2 ~ exp(^-)), is strictly positive at the interval 
(oi, 02) and strictly negative outside of it. Therefore, 
this case doesn't admit the instanton. 

Case 4. n > 1. The potential u(a) is strictly positive. 
The instanton doesn't exist either. 

Both of a newly founded solutions possess only 0(4) 
invariance just like Hawking- Turok instanton (so the Eu- 
clidean region is a deformed four sphere) but, unlike to 
it, they are all non-singular. Note that if the value a is 
sufficiently large then one can neglect the second term in 
(|30p (after the analytic continuation into the Lorentzian 
region) therefore, as in the case of the usual 0(5) instan- 
ton, one can get the de Sitter universe, i.e. the inflation. 

The equation (|2"5|) has no terms with A. In other words, 
the scale factor a(r) doesn't depend on the value A (al- 
though being dependant on the Uq). Therefore, the full 



Euclidean action = S, 



grav 



Sficid has the form, 



S E — So — ASi, 

where So and Si are both independent of the A. Re- 
turning to what has been said in Introduction, there ex- 
ist three common ways to describe the quantum cosmol- 
ogy: the Hartle-Hawking wave function exp(— S E /h) 7 the 
Linde wave function exp(+S E /h) and the tunneling wave 
function. In the second Section we have been working 
with the tunneling wave function. In case of instantons 
situation becomes slightly different. If Si > then (as a 
first, tree semiclassical approximation) we should choose 
the Linde wave function, whereas for the case Si < the 
Hartle-Hawking wave function seems more naturally. 

In conclusion, we note that another choice of C (C < 
and C = 0) eliminates any possible instantons. 
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V. DISCUSSION 

VSL models contain both some of the promising pos- 
itive features [l| and some shortcomings and unusual 
(unphysical?) features as well [T3|. But, as we have 
shown, application of the VSL principle to the quan- 
tum cosmology indeed results in amazing previously un- 
expected observations. The first observation is that 
the semiclassical tunneling nucleation probability in VSL 
quantum cosmology is quit different from the one in 
quantum cosmology with c=const. In the first case this 
probability can be strongly suppressed for large values 
of A whereas in the second case it is strongly suppressed 
for small values of A. This is interesting, although we 
still can't say that VSL quantum cosmology definitely 
results in solution of the A-mystery. The problem here is 
the validity the WKB wave function. And what is more, 
throughout the calculations we have been omitting all 
preexponential factors (or one loop quantum correction) 
which can be essential ones near the turning points. An- 
other troublesome question is the effective potentials in 
VSL models, being unbounded from below at a — ► 0. 
The naive way to solve this problem is to use the Heisen- 
berg uncertainty relation to find those potentials with the 
ground state. However, this is just a crude estimation. 
To describe the quantum nucleation of universe we have 
to find the instanton solution which, being a stationary 
point of the Euclidean action, gives the dominant contri- 
bution to the Euclidean path integral. As we have seen, 
such solutions indeed exist in VSL models. Those instan- 
tons are 0(4) invariant, are non-singular, and provide an 
inflation as well. They describe the quantum nucleation 
of universe from "nothing" and, what is more, upon us- 
age of these solutions we can obtain the probability of a 
nucleation which is suppressed for large value of A (as in 
see Sec. II) using either Linde or Hartlc-Hawking wave 
function. 

Note, that we can weaken the condition n > —1/5 to 
obtain a singular instanton suffering the integrable singu- 
larity (i.e. such that the instanton action will be finite) 
in the way of the Hawking- Turok instanton. However, 
there exist some arguments [14j, that such singularities, 
even being integrable, still lead to serious problems with 
solutions. 

In conclusion, we note that obtained instantons both 
have a free parameter (w 2 ) so we are free to use the an- 
thropic approach to find the most probable values of A 
too. 
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APPENDIX A: THE CASE n > 

Usually those of authors, who apply an idea of vari- 
ability of the light speed to cosmology, are restricting 
themselves to studying the models with negative values 
of n (upon doing so they became able to solve some of 
the cosmological problems without involvement of infla- 
tion). On the other hand, as we have shown in the last 
section, the models with positive n involve a nonsingular 
instanton, and hence are sensible from the point of view 
of quantum cosmology. Moreover, as we shall see, such 
models can also describe our observed universe; in par- 
ticular, they provides a solution to the flatness problem 
and describe the present acceleration of universe. Fi- 
nally, the models with positive n can naturally result in 
Big Rip, as will be shown below. 



1. The dust case: flatness and acceleration 

The integration of ([3]) (with assumption p = wc 2 p) 
results in 



M 



3nks 2 a 2(n - 1 '> 



f a 3(w+i) 47rG(3w; + 2n+l) 4ttG(3w; + 2n + 3) ' 

(Al) 

where M=const. The pressure is 



V 



ws M 3kns wa 



4„„„4n— 2 



s A nkwa An 



4nG{3w + 2n + 1) 4ttG(3w + 2n + 3) ' 



Therefore 

a_ 4ttGM(3w+1) s 2 nk(3w + l)^"" 1 ) 

a ~ 3a 3 ("'+ 1 ) 3w + 2n + 1 

s 2 A(w + l){n+l)a 2n 

+ 3w + 2n + 3 ' 
An equation on the Hubble constant H = a/a is 



(A2) 



2 _ 8irGM (3w + l)ks 2 a 2( - n - 1 > As 2 (w + l)a 2n 

' ~ 3a 3 ( w + 1 ) 3w + 2n + 1 h 3w + 2n + 3 

(A3) 
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If w — then equations (|A2[) and (|A3|1 are reduced to 



4ttGM 
a 



ifca 2 ^- 1 ) s 2 A(n+l)a 2 



if 2 



3a 3 2n + 1 2n + 3 

8ttGM fcsV^" 1 ) As 2 a 2n 



(A4) 



3a 3 2n + 1 2n 

The modern observations result in 



«o 
«o 



7g 2 
10 ' 



(A5) 



where H = h x 0.324 x 10~ 19 s _1 , with h = 45 -f- 75, 
10 28 cm. Thus we have 



a 4?rGAf 
a 



3a 3 



a 2 (l + 2n) 



cgA(n + 1) 
3 + 2n 



and 



Hi 



InGM 
3ag 



(l + 2n)a 2 , 3 + 2n 
Using ![A"5j) . (|X6)) and (|A7| we get to 

3a [agfl^(l + 2n)(3 + lOra) + lOfcc 2 ,] 



M 



407rG(l + 2n)(3 + 2ra) 



A = 



5o 2 r 2 



(A6) 



(A7) 



(A8) 



(A9) 



hence, concluding that such model are indeed able to 
describe our observable universe. 

Let n be positive albeit very small quantity: n < 1 
(and w = 0). In this case the equations ([AT]) . l[A8 |) .(|A4] l 
will be reduced to 



M ns 2 ka 2 



inks 1 



M = 



H 2 = 



12ttG 47rGa 2 ' 



a (3H 2 a 2 + lOfccg) 
40ttG 



3a 3 



a 
a 



AnGM s 2 nk As 2 a 



2 n 1n 



3a 3 



(A10) 



(All) 



It easily follows from the (|A11|) that the evolution of the 
universe can be divided into three stages: 

Stage I. For values of a that are close enough to 
the initial singularity (i.e. a — ► 0) the first term in (|A11|) 
dominates above the second and the third ones. The 
very young universe has been a flat one. [l[ 



Stage II. The second term in (|A11[) dominates above the 
first and third ones. Here we can't consider the universe 
as flat. But we can do it on the last stage. 



Stage III. The last term in (|A11[) dominates above 
the first and second ones. This stage is the one best 
fit to describe the universe we resides in. In order to 
estimate the starting time of this phase we shall use 
(|A11[) . In fact, it is easy to show that last stage will 
start only if 



z < 



-1, 



(A12) 



where z = a^/a— lisa red shift. Therefore, the modern 
astronomical observations can give us the upper bound 
on n. 



1. Big Rip and Big Trip 



The integrals, appearing in (|A10[) (and (jAlip ) in case 
of general position can be extremely hard to solve. How- 
ever, we can make everything a lot more easier, noting 
that at final stage both of the first terms in equation 
(|A10p can actually be omitted. Keeping this in mind, we 
come to 



ait) 




(A13) 



here £*=const. The value of (i.e. a time of a big rip 
occurrence, see ((4])) can be obtained via integration of 
unabridged equation (lAlOj) . Thus, quite surprisingly, we 
get the big rip (See [161 ] about big rip in "usual cosmol- 
ogy" with c =const). 

Remark 3. In the case of "usual cosmology" the 
big rip singularity can exist because of phantom energy. 
There are few ways to escape of future big rip singular- 
ity: (i) to consider phantom energy just as some effective 
models (see [13], [3, 03); (ii) to use quantum effects to 
delay the singularity [20(; (iii) to use new time variable 
such that the big rip singularity will be point at infinity 
(t — > oo ) [ll[; (iv) to avoid big rip via another cosmolog- 
ical "Big": big trip (see below). 

In (22[ Pedro F. Gonzalez-Diaz had shown that phan- 
tom energy can lead to an achronal cosmic future where 
the wormholes become infinite before the arisement of the 
big rip singularity. Soon after that, with the continuing 
accretion of a phantom energy, any wormhole becomes 
the Einstein- Rosen bridge (see also [23|, [13]). Pedro F. 
Gonzalez-Diaz has suggested that such objects can be 
used by an advanced civilization as the means of escape 
from the big rip, but, via usage of Bekenstein Bound we 
have shown it to be impossible due to the very strong re- 
strictions laid on the total amount of information which 
can be sent through this bridge 25]. It is very interest- 
ing that in VSL models with big rip the escape is still 
possible. 
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To show this, we shall use Babichev, Dokuchaev and 
Eroshenko (BDE) equation for the throat radius of a 
Morris-Thorne wormhole [26[ in the dust universe with 
VSL during the last stage of evolution with 

s 2 nAa 2n 



[1] 



12nG 



In this case the BDE equation takes the form 



1 



1 7TW3A 



b(t) b 



log 



t,, — tr 



(A14) 



where b is a throat radius and D is a positive constant. 
Hence, in assumption that the equation of BDE remains 
correct in VSL models, we conclude that there exist a big 
trip at t = T < t», where 



T = U 



(t* - 1 ) exp 



irDbo 



and bo is initial value of b. Now we can calculate the 
horizon distance in such universe: 



R c (t) = a(t) 



c(t')dt' 



it a(t>) ■ 
Substituting (|A13|I into this expression we get 

R c (t) -- 

at the last stage of evolution. Thus, we have no infor- 
mation bound like the one taking place in models with 
phantom fields (where R c ~ * when t — ► i*) and so, the 
advanced civilization can in principle use such objects as 
a means of escape from the big rip singularity. 

Another interesting situation is connected with the 
quantum gravity effects. Such effects take place at a 
scales a ~ Lpi. Since 



^pi — 



Gh Gh {a>o\ 
c 3 Cn V a ) 



3/7 



we can see that Lpi — > oo at a — > (or, in other words, 
the value of red shift z — > oo). It now follows, that the 
universe in vicinity of a Big Bang has been quantum. In 
other words, in order to successfully describe the universe 
for large z one needs the quantum gravity theory. On 
the other hand, close to the big rip we get Lpi — ► 0! 
Therefore, in VSL models, the big rip is a purely classical 
effect and we don't need a quantum gravity (either string 
theory or loop quantum gravity) in order to describe big 
rip in such universe. 
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